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We report a search for the decay B° — > J/ip(j>, using a sample of 657 x 10 6 BB pairs collected 
with the Belle detector at the T(4S) resonance. No statistically significant signal is found and an 
upper limit for the branching fraction is determined to be B(B° — » J/ip4>) < 9-4 X 10~ 7 at 90% 
confidence level. 
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Studies of exclusive B meson decays to charmonium 
play an important role in exploring CP violation [lj and 
establishing the Kobayashi-Maskawa anzatz Q for CP 
violation in the Standard Model. Such studies have also 
resulted in observations of new resonant states that in- 
clude a (cc) pair 

HHH- The deca y B " -> J /^0 

is ex- 
pected to proceed mainly via a Cabibbo-suppressed and 
color-suppressed transition (b — > ccd) with rescattering, 
as shown in Fig. [1] In B decays, effects presumably due 
to rescattering have been seen in various decay processes. 
For example, the large branching fractions observed for 
B° -> DJK+ @ and B~ -> XcqK' Q] decays can be at- 
tributed to rescattering processes @, £ j . An isospin anal- 
ysis on B — » DK^ decays indicates significant final state 
rescattering effects [l(| . Final state rescattering may play 
an important role in understanding patterns of CP asym- 
metries in B decays to two charmless pscudoscalars 11 1 . 
Studies of B decays such as B° — > J/ip<p, which would 
proceed mainly via rescattering, provide useful informa- 
tion for understanding rescattering mechanisms. Previ- 
ously, the BaBar collaboration reported a search for this 
decay mode and set an upper limit for the branching frac- 
tion B < 9.2 x 10 -6 at the 90% confidence level based 
on 56 x 10 6 BB pairs 0]. 
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FIG. 1: Quark-level diagram for B° — > J/tp4> decay 

In this paper, we report the results of a search for the 
decay mode B° — > J/ipcp using the Belle detector [13| at 



on a 605 fb _1 data sample containing 657 x 10 6 BB 
pairs. This sample is more than an order of magnitude 
larger than that used previously. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters (ACC), a barrel- like ar- 
rangement of timc-of-flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL). These detectors are located inside a su- 
perconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside of the 
coil is instrumented to detect K\ mesons and to identify 
muons (KLM). 

Events are required to pass a basic hadronic event 
selection fl5j | . To suppress continuum background 
(e + e~ — * qq, where q = u,d,s,c), we require the ratio 
of the second to zeroth Fox- Wolfram moments [l6| to be 
less than 0.5. 

Candidates for B° — > J/tp<p decays are reconstructed 
from the decays J/ip — ► l + t~ (I = e, /Lt) and (f> — > K + K~. 
The selection criteria for the J/ip decaying to £ + £~ are 
identical to those used in our previous papers [15|. J/ip 
candidates are pairs of oppositely charged tracks that 
originate from a region within 5 cm of the nominal in- 
teraction point (IP) along the beam direction and are 
positively identified as leptons. In order to reduce the 
effect of brcmsstrahlung or final-state radiation, photons 
detected in the ECL within 0.05 radians of the original 
e~ or e + direction are included in the calculation of the 
e + e~{"i) invariant mass. Because of the radiative low- 
mass tail, the J/ip candidates are required to be within an 
asymmetric invariant mass window: —150 (—60) MeV/c 2 
< M e+e -( 7) (M M+M -)-mj /Vj <+36 (+36) MeV/c 2 , where 
mju is the nominal J/ip mass [13] ■ In order to improve 
the J/ip momentum resolution, a vertex and mass con- 
strained fit to the reconstructed J/ip candidates is then 
performed and a loose cut on the vertex fit quality is 
applied. 

In order to identify hadrons, a likelihood Li for each 
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hadron type i (i = ir, K and p) is formed using informa- 
tion from the ACC, the TOF, and dE/dx measurements 
from the CDC. Kaons from the (j> meson are selected with 
the requirement Lk/(Lk + L n ) > 0.7, which has an ef- 
ficiency of 90.0% and a 5.9% probability to misidentify 
a pion as kaon. This requirement is chosen to minimize 
the upper limit expected in the absence of a real signal, 
based on studies of signal and background Monte Carlo 
(MC) events. We reconstruct <j) candidates from pairs of 
K + K~~ candidates, where we require the invariant mass 
to be within ±10 McV/c 2 of the nominal <j) mass fl7j ]. 

B° mesons are reconstructed by combining a J/ip 
with a <p candidate. We identify B° candidates us- 
ing two kinematic variables calculated in the center-of- 
mass system: the beam-energy constrained mass (Mb c = 
V^beam _ ^b) an d tne ener gY difference (AE = E B - 
-Ebcam), where £7beam is the beam energy, and Pb and 
Eb are the reconstructed momentum and energy of the 
B° candidate. We select B candidates within the range 
-0.2 GeV < AE < 0.3 GeV and 5.27 GeV/c 2 < M hc < 
5.29 GeV/c 2 for the final analysis. After all selection re- 
quirements, about 4.9% of the events contain more than 
one B° candidate. For these events, we choose the B 
candidate whose daughter particle cj> mass is closest to 
the nominal value. Finally, a total of 85 candidates are 
selected. 

The dominant background comes from BB events with 
B decays to J/ip- We use a MC sample corresponding 
to 3.86 x 10 10 generic BB decays that includes all known 
B — > J/ipX processes to investigate these backgrounds. 
We find that the dominant backgrounds come from B° — > 
J/ipK*° (892) K~tt + ] and B°/- J/ipK^mO)^ 
K~tt + tt°/~] [18|. In both cases, a pion is misidentified as 
a kaon, and in the latter case, the other pion is missed. 
The former has a peak at AE ~ 0.1 GeV, while the 
latter has a broad peak in the negative AE region. The 
remaining background is due to random combinations 
of J/ip and <f> candidates and does not peak in the AE 
distribution (referred to as combinatorial background). 

The signal yield is extracted by performing an un- 
binned extended maximum-likelihood fit to the AE dis- 
tribution of candidate events. The likelihood function is 



given as 



J2N k 



N 



c 



Nl 



n 



Y^N k xP k (AE l ) 



(1) 



where N is the total number of candidate events, i is the 
identifier of the z-th event, Nj- and P k are the yield and 
probability density function (PDF) of the component k, 
which corresponds to the signal, J/ipKi, J/ipK*°, and 
combinatorial backgrounds. 

The signal PDF is modeled by a sum of two Gaus- 
sians. The background PDFs are two Gaussians for 
the J/tpKi component, a bifurcated Gaussian for the 
J/ipK*° component and a second-order polynomial for 



combinatorial background, respectively. The parame- 
ters of these PDFs are determined from MC simulations. 
We use B° -> J/^K*° decay with K*° -> K~tt + as a 
control data sample to correct for small differences be- 
tween data and MC in the mean and width of the signal 
PDF. The J/ipKi component shape is verified by com- 
paring data and MC events in the K + K~ mass sideband 
region (1.04 - 1.10 GeV/c 2 ), while events in the 5.22 
GeV/c 2 < Afbc < 5.26 GeV/c 2 and K + K~ mass side- 
band region are used to check the combinatorial back- 
ground shape. Possible differences between data and MC 
are included in the systematic errors. 



TABLE I: Summary of the results: upper limits are at the 
90% confidence level. 



Signal yield 


4-6 tii 


Significance 


2.3o- 


Upper limit of signal yield (Ybo) 


9.5 


Detection efficiency (e) 


26.2% 


Upper limit of branching fraction 


< 9.4 x 10" 7 



In the fit, all N k values are free parameters. Figure^ 
shows the AE distribution of the B° — ► J/ip<p candidates 
together with the fit result. We obtain a signal yield of 
4.6 Ij'g events with a statistical significance of 2.3cr. This 
statistical significance is defined as y/— 21n(£o/£max), 
where C max and Co denote the maximum likelihood with 
the fitted signal yield and with the yield fixed to zero, 
respectively. The number of misidentified B° — > J/i/jK*° 
events obtained from the fit is 22.5 _± g and is consistent 
with the expectation obtained from MC simulation incor- 
porating the misidentification probability and the world 
average branching fraction [r?j ]. 

As no significant signal is found for the B° — > J/ip4> 
decay mode, we obtain an upper limit on the yield at 
the 90% confidence level (I90) by a frcquentist method 
using ensembles of pseudo-experiments. For a given sig- 
nal yield, 10000 sets of signal and background events are 
generated according to the PDFs, and fits are performed. 
The confidence level is obtained as the fraction of sam- 
ples that give a fit yield larger than that of data (4.6). 
We account for systematic error by smearing the fit yield 
by the total systematic error described below. We scan 
signal yields and obtain Y 90 = 9.5. 

The corresponding branching fraction upper limit is 
determined with 



B < 



Y 9i 



X Arr 



BB X B(J/l/> 



£+£-) x K+K-) 



Here N BB is the number of BB pairs, and we 



(2) 
use 

the world averages [Uj for the branching fractions of 
B(J/ip -► l + t~) and B{4> -> K+K~). The efficiency 
(e = 26.2%) is determined from a signal MC sample with 
the same selection as used for the data, where a correc- 
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FIG. 2: AE distribution for B° — > J /ip<j) candidates. The curves show the signal (red dashed) and the background components 
(cyan dashed for J/ipKi , magenta dot-dashed for J/ipK*° and green triple-dot dashed for combinatorial) as well as the overall 
fit (blue solid curve). 



tion for muon identification efficiency due to differences 
between data and MC is included (described below). The 
fractions of neutral and charged B mesons produced in 
T(45) decays are assumed to be equal. These results are 
summarized in Table |T] and an upper limit at the 90% 
confidence level is obtained 

B(B° -> J/ip<j>) < 9.4 x ICT 7 . (3) 



TABLE III: Summary of systematic uncertainties on signal 
yield (An) by source. 



Uncertainty Source 


(+cr)An 


(-u)An 


#1(1270) 


1.0 


1.2 


K*° 


< 0.1 


< 0.1 


Combinatorial background 


0.1 


0.2 


Signal 


< 0.1 


0.1 


Total 


1.0 


1.2 



TABLE II: Summary of systematic uncertainties (%) other 
than signal yield extraction (denominator in Eq. 2). 



Uncertainty Source Uncertainty (%) 



Tracking efficiency 


4.2 


Lepton ID efficiency 


4.2 


Kaon ID efficiency 


2.2 


Polarization 


2.6 


J/ip branching fractions 


1.0 


<j> branching fraction 


1.2 


Number of BB 


1.4 


Total 


7.2 



The sources and sizes of systematic uncertainties are 
summarized in Tables [TT] and Mil The dominant sources 
of systematic error in the reconstruction efficiency are 
tracking efficiency and particle identification. Uncertain- 
ties in the tracking efficiency are estimated by linearly 
summing the momentum-dependent single track system- 
atic errors (~ 1% per track). We use control samples 
of J/ip —> £ + £~ and e + e~ — ► e + e~£ + £~ events to esti- 
mate lepton identification efficiency corrections and un- 



certainties. For the J/tp — > mode, we find the 

efficiency for a muon track in the data to be (4.3 ±3.1)% 
lower than that of MC simulation. We correct the effi- 
ciency for this difference and assign a 3.1% uncertainty 
per muon track. For the J/ip — > e + e~ mode, the dif- 
ference between efficiencies in the data and in the MC 
simulation is small, and we assign a 2.7% uncertainty 
per electron track based on their difference and errors. 
We assign an uncertainty of 1.2% per kaon track, which 
is obtained using kincmatically identified kaons in a 
D*+ -> D°n+[D° -> K-TT+] sample. Because of the 
small energy release in cf> — > K + K~ decay, the selection 
efficiency of B° — > J/ip<j> decays depends only weakly on 
the final state polarization. We use an average of the 
efficiencies for fully longitudinally and transversely po- 
larized cases and assign the difference as a systematic 
error (±2.6% including MC statistical error). The sys- 
tematic errors due to signal and background shapes are 
evaluated by varying each of the PDF parameters by its 
uncertainty. We find that the J /ipK\ component uncer- 
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tainty is dominant and that the total systematic error on 
the signal yield is +21.7%/ - 26.1% (Table EHJ. Adding 
all sources in quadrature and conservatively taking the 
larger of the asymmetric errors, the total systematic error 
is estimated to be 27%. As a cross check of the MC effi- 
ciency and analysis procedure, we apply the same anal- 
ysis procedure to the B° — > J/ipK*° control sample and 
obtain B = (1.24 ± 0.01) x 10~ 3 (the error is statisti- 
cal only). This is consistent with the world average 
within its uncertainty and the estimated systematic error 
of the efficiency mentioned above. 

In summary, we have searched for B — > J/4"P de- 
cays. No statistically significant signal is found and an 
upper limit for this decay is determined to be B(B — > 
J/ip<f>) < 9.4 x 10" 7 at the 90% confidence level. This 
result improves upon the previous result [l2T | by about a 
factor of 10 and imposes a more stringent constraint on 
rescattering effects in B° — > J/ipcj) decays. 
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